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ABSTRACT 



Context. Most of the rotational luminosity of a pulsar is carried away by a relativistic magne- 
tised wind in which the matter energy flux is negligible compared to the Poynting flux. However, 
^ I observations of the Crab nebula for instance clearly indicate that most of the Poynting flux is 

^ ' eventually converted into ultra-relativistic particles. The mechanism responsible for transforma- 

^ ' tion of the electro-magnetic energy into the particle energy remains poorly understood. Near the 

equatorial plane of an obliquely rotating pulsar magnetosphere, the magnetic field reverses po- 
larity with the pulsar period, forming a wind with oppositely directed field lines. This structure 
is called a striped wind; dissipation of alternating fields in the striped wind is the object of our 



> 

(N 

00 . study. 



Aims. The aim of this paper is to study the conditions required for magnetic energy release at the 
termination shock of the striped pulsar wind. Magnetic reconnection is considered via analytical 
' methods and ID relativistic PIC simulations. 

Methods. An analytical condition on the upstream parameters for partial and full magnetic recon- 
nection is derived from the conservation laws of energy, momentum and particle number density 
across the relativistic shock. Furthermore, by using a ID relativistic PIC code, we study in de- 
tail the reconnection process at the termination shock for different upstream Lorentz factors and 
magnetisations. 

Results. We found a very simple criterion for dissipation of alternating fields at the termina- 
tion shock, depending on the upstream parameters of the flow, namely, the magnetisation cr, the 
Larmor radius re and the wavelength / of the striped wind. The model depends also on a free 
parameter ^ > 1, which is the ratio of the current sheet width to the particle Larmor radius. It 
is found that for cr » U^r^, all the Poynting flux is converted into particle energy whereas for 
cr <^ (//^ tb)^'"', no dissipation occurs. In the latter case, the shock can be accurately described by 
the ideal MHD shock conditions. Finally, ID relativistic PIC simulations confirm this prediction 
and enable us to fix the free parameter ^ in the analytical model. 

Conclusions. Alternating magnetic fields annihilate easily at relativistic highly magnetised 
shocks. In plerions, our condition for full magnetic dissipation is satisfied at the termination 
shock so that the Poynting flux may be converted into ultra-relativistic particles not in the pulsar 
wind but just at the termination shock. The constraints are more severe for the intra-binary shocks 
in double pulsar systems. Available models explaining observations require low magnetisation in 
the downstream flow. The condition that the magnetic field dissipates at the intra-binary shock 
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implies an upper limit on the pair multiplicity in the pulsar wind k. We found k < few x 10"* for 
PSR 1259-63 and PSR 1957+20. In the double pulsar PSR 0737-3039, the radio emission from 
the pulsar B is modulated with the period of the pulsar A, which implies that the striped structure 
is not erased completely; this gives a lower limit for a- > 3 10. 

Key words. Acceleration of particles - MHD - Shock waves - pulsars: general - Methods: 
analytical - Methods: numerical 



1. INTRODUCTION 



Relativistic shock fronts and currents sheets in relativi stic flows pla y an important role in astro 
physical models of gamma-ray b ursts, for a revi ew see 
and pulsar winds, see for instance 



Michel! (I2005I) and 



PiranI (120051) . jets in active galactic nuclei 



Kirk! (120051) ■ The underlying plasma is prob- 



ably composed of electrons, positrons and/or protons, whose temperature may be relativistic, i.e. 
comparable or larger than their rest mass energy. A shock front is created whenever a fast flow 
encounters the interstellar or intergalactic medium. Relativistic effects become important when the 
post shock temperature is so high that the speed of particles approaches the speed of light c or 
when the bulk velocity of the flow is close to c. Shock acceleration is used to explain the observed 
radiation from gamma-ray bursts or active galactic nuclei. In this paper, we focus on relativistic 
shocks arising form the interaction of the pulsar wind with it s surrounding medi um. For a detailed 
review on theoretical aspects on pulsar wind and plerions, see lLvubarskvl (l2005al) . We briefly recall 
the main issue in this introduction. 

It is widely assumed that most of the rotational energy of a pulsar is carried away in the form of 
an ultra-relativistic magnetised wind. The outflow is dominated by the magnetic field in the sense 
that the energy carried away by the plasma remains small compared to the Poynting flux. This is 
usually described by the magnetization parameter cr, the ratio of the Poynting flux to the particle 
energy flux, which is very large, cr » 1 . Therefore the total power lost by the pulsar may b e 



conveniently estimated by the loss rate of the rotating magnetic dipole in vacuum, (IMichel , 



1991) 



In the general case of an oblique rotator, energy loss can be thought as being shared between the 
steady axisymmetric component and one oscillating with the period o f the pulsar, th e ratio being 
determined by the angle between the rotational and the magnetic axes. lMichell (Il97lb pointed out 
that such waves, which have a phase speed less than that of light, should evolve into regions of cold 
magnetically dominated plasma separated by naiTow hot current sheets. This structure is ca lled a 
striped pulsar wind and has originally been introduced by|( 



Coronitil(ll990l) and 



Michell(ll994l) . 



The structure of this striped wind can be explained as follows. In the aligned rotator, a radially 
outflowing stream of particles opens up the dipolar magnetic field Unes crossing the light cylin- 
der. Asymptoticall y, the stream fi nes become radial and can be described by the so-called split 
monopole solution (lMichellll973h . It consists of two half magnetic monopoles with equal magni- 
tude but opposite sign separated by the equatorial plane. When crossing this equatorial plane, the 
polarity of the magnetic field reverses, implying a surface current density, i.e. the current sheet. If 
the rotation axis is tilted with respect to the magnetic moment, the current sheet oscillates around 
the equatorial plane and develops a wave structure, expanding radially at a speed close to the light 
velocity. In a radial outflow, the amplitude of the oscillations grows linearly with radius and at large 
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distances can be approximated locally by spherical current sheets separating the stripes of magne 
tised plasma with opposite magnetic polarity. The exact solution for t he oblique split monopole in 



the ide al MHD approximation has been found by iBogovalovl (Il999b . Simulations by ISpitkovskv 



( l2006h in the force-free limit show that beyond the light cylinder, the wind from an oblique dipole 
pulsar magnetosphere is similar to that from the split monopole. The structure of the current sheet 
in the striped wind is shown in Fig. [1] Note that the oscillating current sheet is found in the solar 
wind so that this structure is generic. 




Fig. 1. Structure of the current sheet in the oblique rotator. The rotating neutron star is located at the 
centre. When crossing the current sheet, the magnetic polarity is reversed. Moreover, this structure 
is expanding radially outwards with a constant speed close to the light velocity. 



While expanding into the nebula, the pulsar wind terminates at a standing shock located at a dis- 
tance where the confini ng pressure of the nebula balances the rarn pressu re of the wind, this is called 



the termination shock, dRees & Gunn , 



1974 



Kennel & Coroniti , 



1984 



I; 



Emmering & Chevalier , 



19871). At the shock front, the energy in the wind is released into ultra-relativistic particles respon- 
sible for the observed radiation. Observations of the interaction between the wind and the neigh- 
boring nebula clearly indicate that the electromagnetic energy of the wind is largely converted into 
particle kinetic energy. However, the conversion mechanism and the corresponding acceleration 
process still remain unclear. 



Indeed, the analysis done by different authors. (ITomimatsu , 



Chiueh et al 



1998 



Bogovalov & Tsinganos , 



1999; 



1994 



Lyubarsky & Eichlei 



Beskin et al 



1998 



2001 ) shows that in an 



ultra-relativistic, radial wind, the electric force compensates the magnetic tension so that the flow 
could not be accelerated significantly. This means that the electromagnetic energy is not transferred 
to the plasma. To sum up, in the region where the wind is launched, the magnetisation should be 
high, cr » 1, whereas just beyond the termination shock, it should be small, cr <s; 1. But according 
to the previous results on MHD wind collimation and acceleration, a significant decrease in cr is 
forbidden. Therefore, there is a contradiction which is known as the cr-problem. A promising alter- 
native solution to convert the Poynting flux to particle kinetic energy is investigated in this paper. 
We demonstrate that, under certain conditions, magnetic reconnection in the striped wind occurs 
when it crosses the termination shock. 



Both obser vations of the inne r region of the Crab nebula done by IWeisskopf et al 



(l2000t) and 



the solution by iBogovalovl d 1 9991) show that most of the energy carried by the wind is transported 
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in the equatorial plane of the pulsar wind. In this region, energy is transferred predominantly by the 
alternating magnetic field. So dissipation of alternating fields i n the striped win d could be the main 
energy conversion mechanism in pulsars. It was recognised by 



Michell (1 19821) that 



Usovl(ll975h and|] 

the amplitude of the magnetic oscillations decreases with the distance as r ' whereas the particle 
number density sustaining these waves falls off as where r is the radius in spherical coordinates. 
At a given point, the charg e carriers becom e insufficient to maintain the required current and the 
alternating field annihilates. ICoronitil (1 1 9901) considered magnetic reconnection in the striped wind 
and was lead to the same conclusion. However, the flow is accelerated during the process of mag- 
netic dissipation, which dilates the timescale of the wave d ecay such that the magnetisation remains 
high at the termination shock, dLvubarskv & Kirkl l200lh . Therefore, the wind enters the termina- 



tion shock still dominated b y Poynting flux unless the annihilation rate is nearly to the causal limit, 
dKirk & Skiaeraasenll2003h . However, when the flow enters the shock, the plasma is compressed, 
which could result in the forced annihilation of the alternating magnetic fields. In this paper, we 
address annihilation of alternating fields at a relativistic shock in electron-positron plasma. 



In iLvubarskvl (l2003h . the reconnection of the magnetic field at the termination shock was 
studied phenomenologically by introducing a fraction tj of the magnetic energy dissipated 
at the shock. Relativistic MH D shock fronts with dissipation have aheady been studied by 
Levinson & van PuttenI (1 1997b . The jump conditions for a relativistic perpendicular ideal MHD 
shock with arbitrary magnetisation has been investigated bv lZhang & Kobavashil (l2005h in the con- 
text of GRBs. In the present paper, we use a kinetic description of the plasma in the striped wind. 
First we find jump conditions assuming that the thickness of the current sheets downstream of the 
shock is scaled as the particle Larmor radius; this gives us the fraction of the dissipated magnetic 
energy as a function of upstream parameters and of a phenomenological parameter, ^ > 1, which is 
defined as the ratio of the sheet width to the Larmor radius. Then we perform particle-i n-cell (PIC) 



Lvubarskv 



simulations of the shock in the striped wind (such simulations were shortly described in 
(l2005bl) ) and find the parameter ^ from the simulation results. We show that the alternating fields 
are easily annihilated at the shock front so that the electromagnetic energy of the pulsar wind could 
be readily converted into the plasma energy at the termination shock. This could provide a solution 
to the cr-problem. 

The paper is organised as follows. In Sect.|2] we present the full system for the jump conditions 
of the averaged quantities in the striped wind. In Sect[3] we solve the system analytically. We derive 
an analytical condition for magnetic annihilation for an ultra-relativistic and strongly magnetised 
flow. The results are then checked and extended by numerically solving the jump conditions at the 
termination shock for the average parameters. In SectH) we perform several ID PIC simulations of 
the striped wind by varying the parameters in the upstream plasma hke the Lorentz factor and the 
magnetisation. Two typical situations are presented, the first one demonstrating that no magnetic 
energy is released at all and the second one showing full dissipation of the Poynting flux which is 
converted into particle heating. An empirical law for dissipation is then derived from the full set 
of PIC simulations. The results are compared to those obtained in Sect. [3] In Sect|5] the results are 
applied to pulsars in binary systems. On one hand, magnetic dissipation at the termination shock 
implies an upper limit on the pair multiplicity factor k. Applications are presented for two binary 
pulsars, PSR 1259-63 and PSR 1957+20. On the other hand, an absence of significant dissipation 
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imposes a lower limit on k. This is applied to the double pulsar PSR 0737-3039. The conclusion 
are presented in Sect. |6l 

2. JUMP CONDITIONS IN THE SHOCK 

In this section, we derive the criterion for dissipation of alternating field at the shock by considering 
the jump conditions for the averaged quantities in the flow, namely, the conservation of the energy, 
momentum and the particle number density. We assume that the striped structure survives in the 
downstream flow and find the condition of the total dissipation of the alternating fields from the 
condition that the thickness of the current sheet downstream of the shock approaches the thickness 
of the stripes. 



2.1. Description of tine wind fiow 

Let us consider a one dimensional striped pulsar wind entering the termination shock. By conven- 
tion, we refer to quantities in the shock frame with unprimed letters, while in the proper frame 
of each plasma, quantities are denoted by a prime. However, thermodynamical quantities such as 
pressure p, temperature T, internal energy e and enthalpy w, are always expressed in the proper 
frame, so dropping primes should not lead to any confusion for these quantities. 

The striped wind propagates in the x-direction with relativistic velocity and possesses an alter- 
nating magnetic field directed along the z-axis. Moreover, in this and the next sections we assume 
that its average over one period of the striped wind vanishes. Quantities upstream, i.e. before cross- 
ing the shock discontinuity are subscripted by 1 whereas quantities downstream are subscripted 
by 2. 

Let us describe the state of the incoming plasma. In the shock frame, the upstream Lorentz 
factor of the ultra-relativistic wind is Fi » I. The current sheets are made of a hot, unmagnetised 
plasma with density Whi and a temperature T^i. The distance separating the middle of two succes- 
sive current sheets is denoted by li (half a wavelength of the wind) and their thickness Ai is much 
less than half a wavelength, Ai <k li. The magnetised part of the wind is cold, Td = 0, has a 
density n^i and a magnetic field strength Bi. Moreover, the wind is strongly magnetised such that 
the magnetisation parameter upstream defined by 

" ~~P2 = ~F 2 

/^orfwci ^l()Tln^xmc^ 

is very high, cri » L The enthalpy is simply given by the rest mass energy of the cold ultra- 
relativistic particles, Wci = n^j mc^. The speed of light is c and m is the mass of the leptons 
composing the wind, actually electrons and positrons. We neglect the enthalpy contribution from 
the hot current sheets because their thickness is assumed to be very small. 

Downstream, the wind is decelerated to a Lorentz factor r2 and compressed such that the dis- 
tance separating two successive current sheets shrinks to a length I2. When compressed, the cold 
magnetised component of the wind is heated to a temperat ure T given in the Unfit of high mag- 
netisation, cTi » 1, and high Lorentz factor, Fi » 1, by dKennel & Coronitil Il984al) (see also 
appendixlA]) : 

^.lll (2) 

mc^ 8 F2 
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ks is the Boltzmann constant. We assume here that the shock width is much less than the wavelength 
of the wind. When the shock is between the sheets, ideal MHD shock applies and Eq. (|2]) is locally 
valid. Note also that the shock velocity in the cold part is close to the speed of light c and different 
from the shock velocity in the hot part. Quantities are evaluated in the frame where the shock is at 
rest on average. Eq. ^ is vaUd in any frame. 

Note that even if the upstream flow is a pure entropy wave with a constant magnetic field 
between the sheets, fast magnetosonic waves should be generated beyond the shock. Therefore the 
structure downstream is not steady in the proper frame; there should be oscillations. We consider 
quantities averaged over the wave period and neglect contribution of these oscillations into the 
fluxes. 

Pressure balance between gaseous and magnetic part does therefore apply on both sides of the 
discontinuity. The gaseous pressures are given by : 

Pi = "hi^BThi (3) 

P2 = nl^2^BTh2 (4) 



8 r2 

Note that the "cold" part of the wind is heated to relativistic temperature downstream and therefore 
also contributes to the gaseous pressure via the term p^a^ where the temperature is given according 
to Eq. (|2|l. In other words, upstream we have 

2/uorf 

whereas downstream, taking into account the heated cold component, we find 

Pi = Pel + (7) 

Due to the Lorentz length contraction, the densities in the proper frame and in the shock frame are 
related by n = T «'. On the other hand, the Lorentz transformation of the magnetic field is B = T B' . 



2.2. Jump conditions for average quantities 



We now write down the MHD jump conditions for the qua r itities averaged over one wavelength of 
the wind. Following the procedure described in iLvubarskvl (l2003b for a perpendicular MHD shock, 
the conservation laws can be cast in a form similar to those for a relativistic hydrodynamical flow. 
Noting that in the cold magnetised part, the field is frozen into the plasma 

B\ B2 Bi^ _ Bc2 



(8) 



Fi n^j r2 "<:l "c2 

such that the magnetic field becomes a function of the density, one can introduce the eff^ective 
pressure and enthalpy by : 



Pi 



Pi + 



Wi + 



I 



(9) 
(10) 



Now the jump conditions at the shock discontinuity are given by 
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- conservation of particle number density : 

<ri Pi n\ ) = <r2/32 n2> = </8i m) = </62 n2>; (11) 

- conservation of total energy : 

<r2^inVi) = <r2/32W2); (i2) 

- conservation of total momentum : 

{Y\p\^^'^+v{)^{Y\p\^^'^+v^). (i3) 

Here j0i,2C are the upstream and downstream velocities, correspondingly and the angular brackets, 
(), mean averaging over the wave period, which means spatial averages. We assume that the down- 
stream flow is settled into an equilibrium pattern moving with constant velocity, i.e. that only an 
entropy wave is presented downstream. Then time averages are expressed via spatial averages; e.g. 
the average particle flux is 

1 If' y6 r' 

— I nfidt—— I ndx—— I ndx. (14) 
T Jo T Jq I Jo 

As the time-averaged particle flux (the left-hand side of Eq. (fT4l l) is constant and the total number 
of particles within the wavelength ( Jj^' n dx) is conserved, Eq. ( fT4l i yields a useful relation between 
the upstream and downstream wavelengths: 

f = f (.5, 

Of course, our assumption is not exactly correct; when an entropy wave impinges on the shock, 
both entropy and fast magnetosonic waves generally appear in the downstream flow. One could 
assume that the shock width is small as compared to the sheet width; then one could apply Rankine- 
Hugoniot relations locally and find the time-dependent structure of the downstream flow. However, 
such an assumption is too stringent; it could be valid in fact only when dissipation is negligible. 
In this paper, we look for a criterion for complete dissipation of the alternating magnetic field. 
This criterion is found by extrapolating our analytical model to the case when the downstream 
Larmor radius becomes comparable to the wave period so that the structure assumed in the model 
is already destroyed. Therefore this criterion is anyway very rough. We believe that such a rough 
criterion should be independent of the fine structure of the downstream flow, like fast magnetosonic 
oscillations. PIC simulations indeed show that our criterion for dissipation is roughly correct so that 
our model is viable. Detailed study of the shock interaction with the current sheets in case of partial 
dissipation, as well as generation of fast magnetosonic waves in the downstream flow, is the subject 
of our future work. 

In the hot phases of the plasma, we use the ultra-relativistic equation of state such that w = 4p. 
Therefore, the average effective enthalpies upstream and downstream are 



<nV2) = 4p262 + (l-62) 



, 2 Bi 
n^, mc + 



B 



2 
2 



(16) 



po Ti 



(17) 

2- 

For the remaining of this paper, we introduce the relative thickness of a current sheet by 5, - Ai/li, 
(i - 1,2). Making use of the pressure balance condition Eq.®, one can express the upstream 
enthalpy via the magnetisation parameter Eq. ([!) as 

<^i> = n' mc^ [l-6i+(l+ 6i)cri]. (18) 
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The downstream pressure balance Eq.(|7]i could be also expressed via cri taking into account that 
between the current sheets, the magnetic field is frozen into the plasma, Eq. (jS): 



P2 = g p^'\.2"^C (1+4 — 0-1 



rici 



(19) 



(20) 



Then the downstream enthalpy is written as 

2 r2 

Now the conservation of particle number density, energy and momentum simplifies respectively 
into : 



<^2) = -^<2"^c2 



1 +(1 +^2)2 — CTi 



ySi [(1 -(Ji)nci +<Jinhi] 
2/3in,i [1 -(5, +(1 +(5i)cri] 

IV 1 



Pi [(1 - (52) nc2 + ^2 nh2] 



y82n. 



c2 



1 +(1 +(52)2 — (Ti 



= «c2 ir2y82 



1 +(1 +^2)2 — crj 



1 nc2 CTi 

H 1 

4r2 ncir2 



(21) 
(22) 

i£i 

2 r2 

(23) 



The striped wind downstream is described by four parameters, namely, the particle number density 
in the hot (unmagnetised) and cold (magnetised) part, n^a and «c2, the Lorentz factor r2 and the 
relative current sheet thickness 52- The available three equations ( 1211 1. ( l22l i. ( |23] ) should be comple- 
mented by one more equation; some assumption about microphysics of the reconnection process is 
necessary to close the system. 

The most natural assumption is that in the course of the reconnection process the current sheet 
thickness scales as the Larmor radius of the particle in the sheet; then the system is closed by 
introducing a parameter ^ > 1 such that the downstream current sheet thickness in the proper frame 
is defined by 



A2=^4 

The Larmor radius in the downstream plasma frame is given by ; 



(24) 



(25) 



« \q\B'^c \q\B2C 

q is the charge of a lepton (electron or positron), q - +e. Expressing the temperature via the 
pressure Eq. ^ and substituting the pressure balance condition Eq. ( fT9] l. one can write Eq. (l24l) as 



a; = 5- 1 1 + 4 — 0-1 — ro; 

"cl / «h2 



(26) 



where we introduced the Larmor radius related to the upstream bulk velocity Px c and defined by 
T\mc 



\q\Bi 



This quantity could be deduced from the upstream parameters as 



1 

'"0 = 7-: 



Fi m 



\q\ \t 



(27) 



(28) 



flO CTl Hcl 

Note that the expression (|27] | is not the true Larmor radius of the particles in the upstream flow, the 
last being dependent on the plasma temperature. Transforming to the shock frame, A2 - A'2lT2, 
dividing by h and making use of Eq. ( fTsT i, one can finally write the closure condition as 



02 - - I 1 + 4 O"! 

"el / nh2P2 ^ih 



(29) 
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3. SOLUTIONS OF THE JUMP CONDITIONS 



9 



The set of equations (12111 -(l22ll-([23Tl-(|29]l is first solved analytically under the assumption that the 
fraction of the dissipated magnetic energy is small so that the downstream current sheet remains 
narrow and the downstream flow remains ultra-relativistic and strongly magnetised. Extrapolating 
of the obtained asymptotics to the case when the thickness of the current sheet in the downstream 
flow becomes comparable to the wavelength, we obtain an analytical criterion for the full dis- 
sipation of the alternating magnetic fields at the shock front. More general conditions are then 
investigated by solving the system numerically. 

3. 1 . Analytical asymptotic solution 

We are interested in a highly relativistic, Poynting dominated upstream flow so that crj » 1, 
Fi » 1,^1 <K 1 . I n the ideal MHD, such a flo w would remain highly relativistic downstream of the 



shock, r2 - -\/cr7. lKennel & Coronitil(ll984ah and appendixlAl When some fraction of the magnetic 



energy is dissipated at the shock front, the velocity of the d ownstream flow de creases and reaches 



20031). In this subsection. 



^2 - 1/3 when the magnetic field is completely dissipated (ILyubarskyl 
we assume that only a small fraction of the magnetic energy is dissipated so that downstream of 
the shock, the current sheets remain narrow, ^2 ^ 1 and the flow remains ultra-relativistic, r2 » 1 . 
We solve Eqs. (l2TI)-(l22b-(l23l)-(|29t assuming that 1 /cri, l/Ti ^ and 5i_2 are small. In the zeroth order 
approximation in these parameters (when all of them are equal to zero), the three equations (1211 1- 
(l22] )-(|23l) are reduced to the same equality 

tici - nc2 (30) 

so the system is nearly degenerate. It is therefore necessary to retain higher order terms. For a while 
we do not assume any relations between these parameters and just expand Eqs. (l2Tl i- (l22l i-(l23Tl to 
the first non-vanishing order in each of them independently. In the obtained equations, both the 
zero and the first order terms are presented. Dealing with equations containing terms of different 
order is difficult. One can simplify the problem if one uses Eq. ( l22b in order to eliminate the zeroth 
order terms from Eq. (|2TI ) and Eq. ( |23l l. correspondingly; then one gets two equations containing 
only small order terms; these equations could be complemented by the zeroth order equation ( |30] | 
and also by Eq. ( |29l ) in the zeroth order. 

Expanding Eq. ( |2T] ) to first non-vanishing order in small parameters we find 

Wcl Wc2 

"ci - nc2 = Si (Wci - "hi) - S2 (nc2 - nh2) + ^^2 ~ 

We put the zeroth order terms to the left hand side and the small terms to the right hand side; one 
can see that the difference n^i - n^2 is small. Now let us divide Eq. (l22l l by cr\ so that the leading 
order terms be of zeroth order; then expanding in the small parameters to the first non-vanishing 
order yields 
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One sees that neglecting small order terms would result again to Eq. ( [30l l therefore the small terms 
should be retained. Substituting Eq. (ISTT l into the zeroth order terms and Eq. ( l30t into the rest of 
the terms, one gets 
111 2 2 

Here, we introduced the fraction of cold to hot particle densities by 

Z = — (34) 

"A 

In order to eliminate the zeroth order terms from Eqs. (|22]) and (l23T l. one can simply extract 
one of them from another because the zeroth order terms in these equations are the same (those 
containing Y^crx). Making use of Eq. ( |30] | in the rest of the terms, we arrive at 



\ ^ _ £i 



1 1 



(35) 



The setofEq. (I30ll-(l33ll-(l35ll is equivalent to Eq. (|2B-(l22li-(|23]l if the parameters 5i, l/cri and l/E^ 
are small. 

Recall that the speed of the shocked plasma cannot exceed the speed of the fast magnetosonic 
wave having a Lorentz factor Ffms - -^cr\j{\ - c\) where Cs is the sound speed. For an ultra- 
relativistic gas, Cj ~ 1 / V3 and therefore r2 < Ffms ~ V3 o"i/2. Thus, for a super-magnetosonic 
upstream flow satisfying Fi » -y^' downstream Lorentz factor always satisfies F2 <s; Fi. 
Assuming also that upstream of the shock, contribution of the hot plasma in the sheet to the total 
plasma energy is small such that 5\ <K 1/cri, one reduces Eqs. ( |33] | and ( |35] | to 



1 1 2 

&i (3 - — ); (36) 

2 



2o-i 2F? Z2 



+ = — T- (37) 

4o-i 4F? 



1 

A 1 

- 2 

Combining Eqs. (|37] | and (|36] | one immediately gets Z2 =2/5. 

The two remained unknowns, ^2 and F2, could be found from Eq. (ITTT l and the closure condition 
Eq. ( |29] l. In the limit cri » 1, Fi 2 » 1 the last is written, with account of Z2 - 2/5, as 

,2=^^^. (38) 
J 1 2 'I 

Making use of Eq. (|37] | and Eq. (l38T l. we can now analyze dissipation of the alternating magnetic 
fields at the shock front. 

An interesting quantity is the ratio of the downstream magnetic energy flux to the downstream 
matter energy flux. This is the true magnetisation parameter of the shocked flow defined as 

(1 - '^2)g^ 2(l-g2)«c2Cri 

0"2 — r — (39) 

4/^o[52;'2 + (l -5i)Pcl¥\ «cl +4^2nc2 0"l 

If the dissipated fraction of the magnetic field is sufficiently small, the flow should satisfy the ideal 
MHD jump conditions; then the downstream magnetisation is given by cr^ - Icrx ( appendix lAb. 
Inspecting Eq. ( [39] ), this condition leads to 

(52 « (40) 
4cri 

Eq. (l40t means that contribution of the plasma in the current sheets to the plasma energy in the 
downstream flow remains negligible small so that magnetic dissipation does not affect the dynamics 
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of the flow. Moreover at the condition Eq. ( |40] i. Eq. dJTj i yields r2 = which is the usual result 
for ideal perpendicular MHD shock. Substituting this result to Eq. dSST l gives 



„ . 1 Vo^n) 



(41) 



Now the condition ( |40] l for the negligible dissipation can finally be written in terms of the upstream 
parameters as 



(Ti <sc 



5/i 
4^ro 



2/3 



(42) 



Now let a significant fraction of magnetic energy dissipate. Then the condition opposite to 
Eq. ( |40] | is fulfilled so that one can neglect the second term in the left hand side of Eq. (|37] |. 
Together with Eq. ( [38] ) this yields 

jSi 4^ a\ ro' 

\2 



(43) 



62 = 



2 ^ 0-1 rp 
5h 



(44) 



At the condition opposite to Eq. (|40| ). Eq. ( |43] | gives r2 < yoT- In order to have physically mean- 
ingful results, we assume that ^2 < 1, i e- that the current sheet thickness can never exceed the 
distance between two stripes. Full dissipation occurs when the width of the current sheets becomes 
as large as the distance between two successive sheets, which means (^2 — ^ 1 ■ This occurs when 
the expression in the curly brackets in Eq. (l44l l goes to unity. Now the condition of full magnetic 
dissipation could be written as 
5/1 



0-1 > 



2^ro 



(45) 



When this limit is reached, the Lorentz factor downstream, Eq. ( |43] ). becomes non-relativistic 
r2 ~ 1. More precisely, setting ^2 = 1 in the exact jump conditions Eq. (l22b and ( |23] ). we get 



2Pi n,i (1 +0-1) = j62 nc2 1 1 + 4 — cri 

rici 



j = nc2 |r2y62 


1 ^ "^2 
1+4 (Ti 







+ 



1 



ricio-i 



4r2 HciTz 



in 

2 r2 



In the highly magnetised and ultrarelativistic hmit, we find 
1 



/32 + 



4r2^2 



= 1 



(46) 
(47) 

(48) 



The only physically acceptable solution is ^62 = 1/3. Therefore, in case of full dissipation, the 
downstream parameters are the same as in the non-magnetized shock. 
We summarise the aforementioned results in Table [T] 





0-, » 5/i/(^ro) 


1/5 « /i /(fro 0-1) « Vo=r 


0-1 «[5h/(4^ro)f" 


Reconnection level 
S2 

T2 


Full 
ai 1 

« 3/(2 V2) 


Partial 

(2ffriro/5/,)2 
5/i/(4fcriro) 


Negligible 
« 1/4 0-1 



Table 1. Summary of the analytical criterion for the dissipation of alternating fields at the shock. 
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3.2. Numerical solution to the jump conditions 

In order to investigate the reconnection properties in less restricted limits than those used to obtain 
an analytical solution to the jump condition in Table [T] we solve numerically the average conser- 
vation laws of particles, energy and momentum, Eq.(|2TJ, (l22l l. (l23T l supplemented with Eq. (|29] l 
for the current sheet thickness. We remove the assumption cti » 1 but keep a highly relativistic 
supermagnetosonic flow, Fi » -^(r7, with thin current sheets, 5i <k l/cri. 

We solve the full system of Eq.(l2Tli. (|22] |. ( l23T l. Q and ( l26b with upstream flow conditions 
given by Ti = 10"*, rid = "hi = 1, 5i = 10"** and different values for the parameter ^ = 10, 10^, 10"'. 
Note that if the flow is super-magnetosonic, Fi » yfa^, and the contribution of the hot plasma in 
the sheets to the overall plasma energy is negligible. Si <K l/tri, the results are independent of ^i, 
Zi and Fi; the only important quantity being the upstream magnetisation cri. The free parameter^ 
will be determined later by performing PIC simulations, see Sect. |4] 

Results for different upstream magnetisation cri and diff'erent half wavelength of the striped 
wind 1 1 are summarised in Fig.|2] The results scale with ^ as expected so we divide the abscissa by 
the parameter ^ to get "universal" curves. When the contribution of the hot plasma in the current 
sheet to the overall plasma energy downstream is small, i.e. for ^2 0"! <K 1, dissipation is negligible 
and the ideal MHD jump conditions apply. Inspecting Fig.|2^), we conclude that this should happen 
whenever tri < (h/^ ro)^^^ so that we rec over the no dissipation crit erion Eq. (l42T i. Therefore, the 
Lorentz factor is equal to F^ x cri +9/8 (IKennel & Coronitil (Il984ah . appendixlAll as readily seen 
in Fig. I2J)). Because the upstream flow is not necessarily highly magnetised, in order to achieve 
satisfactory accur acy, we kept two terms in the expansion of the Lorentz factor F2 with respect 



to cTj as given in 



Kennel & Coronitil ( Il984ah . In this region, the density of the cold part is close 
to the second order approximation «c2/«ci ~ (1 + 1/2 cri - 3/16crj). The wind remains entirely 
dominated by the electromagnetic energy flux, 0-2 ~ 2 cri. The dissipation is complete when 62 
approaches unity. Fig. |2};), so that the striped structure is removed. Equivalently, it corresponds 
to the case where <J2 <K 1. Inspecting the Fig. |2f), we expect this to happen when cri > li/^ro; 
so we retrieve the condition Eq. (05]) for full Poynting flux dissipation. In this case, the magnetic 
field is converted into particle heating. The thermaUsation is complete and the flow downstream is 
purely hydrodynamical. The total energy flux is entirely carried by the matter energy flux, 0-2 «: 1 . 
Note that the region where (^2 > lis meaningless. This can imply some negative magnetisation for 
low (Ti as clearly seen in Fig.|2^). The fraction of cold to hot particle densities Z2 is close to 2/5, 
fig|2]a). Actually, for sufficiently large magnetisation cri, it does not vary much and remains close 
to this value, let reconnection be or not. The Lorentz factor becomes close to unity F2 ~ 1, fig|2]b). 

To resume, the numerical solution of the jump conditions confirm our analytical expectations 
summarised in Table [T] Results have been extended to less restricted upstream flows, the magneti- 
sation being not necessarily very high. A good estimate for Z2 is 2/5. It is a good approximation 
whenever cri > 100, as depicted by a solid thick line in Fig.|2^). 

In Fig. |2}l), we plot the fraction of cold to hot plasma component given by 
1 -62 

Zi (49) 

02 

It is clearly seen that this ratio is not constant and much less than the upstream fraction given by 
Zi = 10^ (50) 
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Fig. 2. Several downstream parameters obtained by numerical solution of the jump conditions. The 
ratio Z2 = rtci/rihi is shown in Fig a), the downstream Lorentz factor in Fig b), the downstream 
relative current sheet thickness in Fig c). The fraction of cold to hot plasma component after shock 
crossing, evaluated by Eq. ( |49l ) is shown in Fig d). Even in the low dissipation limit, it is not equal 
to the upstream fraction, denoted by solid thick line at 10^, see Eq. ( fSOb . The ratio of downstream 
to upstream magnetisation 0-2! o-\ is shown in Fig e), the downstream magnetisation in Fig f). Note 
that for the ratio 0-2! cr\ , the abscissa scales like cr"^^^ in order to overlap all curves independently of 
o"! and to match condition ( l42l l. All quantities except for fig e) are plotted versus the adimensional 
parameter Zi crj ro. Different colors associated with different symbols represent different logjQ ctx 
as detailed in the legend. The asymptotes in the limit of no dissipation (ideal MHD shock) in fig. b) 
are shown in colored dotted lines. 

and shown by a solid thick line. The reason is that some magnetic dissipation occurs even at the 
condition (l42T l. The amount of the energy dissipated is not sufficient to affect the jump conditions 
but some particles from the cold part diffuse across the magnetic field lines and enter the sheets. 

In a final step, we determine the free parameter ^ with help of PIC simulations. This is discussed 
in the next section. 



4. TERMINATION SHOCK: PIC SIMULATIONS 



We desig ned a fully relativistic and electromagnetic ID PIC code following the algorithms de- 
scribed in lBirdsall & Langdonl(l2005b . Particle trajectories are advanced in time by integrating the 



relativistic equation of motion due to the Lorentz force. The longitudinal electric field is found by 
solving Poisson equation whereas the transverse component of the electromagnetic field are com- 
puted by the remaining component of Maxwell equations introducing a left and right-going wave 



14 Petri and L5aibarsky: Magnetic reconnection in a striped pulsar wind 

such as F"^ - Ey+c B^, see lBirdsall & Langdonl(l2005b . The simulation is one dimensional in space 



along the x-axis and two dimensional in velocity in the plane (xOy). 

The striped wind propagates along the jc-axis and hits a wall located at the right boundary of the 
simulation box of length L. We also impose no incoming electromagnetic wave at this boundary. 
The particle momentum vector possesses a longitudinal component and a transversal compo- 
nent py. They propagate from the left to the right and hit the solid wall at x - L. This means that 
particles are reflected at this boundary, the p^ momentum is reversed whereas the py component 
remains unchanged. The magnetic field is directed along the z-axis and reverses polarity when 
crossing each current sheet. Electromagnetic waves leave the simulation bo x at the right bound ary 



without any reflection. Similar simulations have already been performed bv lLvubarskvl ( I2005H) . 

A typical initial situation showing the magnetic field structure, the lepton distribution functions, 
the temperature and the drift speed is shown in Fig. [3] Close to the right wall, we impose a decreas- 
ing magnetic field strength B^ix) which almost vanishes at x - L, Fig. |3^). The magnetic polarity 
reversal is smooth, it is not squared but evolves according to the following expression 

B/x) = Bii tanh[ A (a + b co&(A x))] D(x) (5 1 ) 

where Bo is the maximal intensity of the magnetic field, a, b, A and A are constants, A prescribing 
the thickness of a current sheet and A the number of periods in the simulation box. T his particular 
expres sion is dictated by the structure of the striped wind, see for instance Eq. (1) in 



Petri & Kirk 



(|2005|). D(x) is a function introduced to decrease the magnetic field intensity close to the wall, at 
X - L. Moreover, because magnetic pressure is balanced by gaseous pressure, the temperature T(x) 
in the gas also decreases close to the right wall where it almost vanishes. Fig. |3}l), in accordance 
with the magnetic field behavior Lower temperature implies weaker spread in particle momentum 
space, Fig.ISj)) and c) (note that we keep track of only 50.000 particles, chosen randomly at each 
time step, in order to avoid too large data files). We also take the bulk Lorentz factor of the flow, T\ , 
decreasing near the right wall. The average longitudinal momentum p^ also decreases in relation 
with the decreasing Fi . With this initial condition, we avoid the formation of too large a transient 
when the plasma first collides with the right wall. 

The plasma moves in the positive x direction with the bulk Lorentz factor Fi - 20. The temper- 
ature in the plasma is obtained from the equilibrium conditions, namely pressure balance between 
gaseous and magnetic part. In the current sheet, the temperature is much higher than in the mag- 
netised part, Fig.[3}i). This implies a much larger spread in momentum as seen in Fig.lSj)) and c). 
The average magnetic field a -< B > /Bq is not necessarily zero because two successive stripes 
generally have diff'erent width but the same magnetic field intensity Bo- In the pulsar wind, the av- 
erage magnetic field only vanishes in the equatorial plane, 6 = n/lin spherical coordinates (r, 6, ip). 
Actually, in Fig. [3^) as well as in the simulation results presented below (Fig.|4^) and Fig.[8t)), the 
average magnetic field is o- » 0.12. We choose an average value different from zero not because 
of numerical stability requirement but because we want to demonstrate that even in case of full 
dissipation the average magnetic flux downstream the shock is preserved. We also performed a set 
of simulations with zero average magnetic field. The results obtained are not aff'ected and remain 
qualitatively the same in both cases as wiU be shown later Moreover, passing from one current 
sheet to the next one, the drift speed for each species reverses sign. Fig. [3^). Indeed, in a current 
sheet, the magnetic field, oriented in the z direction, is sustained by an electric current flowing in the 
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y direction. This current is generated by electrons moving in, let's say, the positive y direction at a 
speed Us whereas positrons are moving in the opposite direction at a speed -Us- Charge neutrality 
is maintained but a net total current exists. 




2 .J 





.1^ 



— i7i — »-r-^:i ' - i 



■ JuJIa. 



ML 



-.'1. 



□ D.I v.ia B,* m,m%' 



Fig. 3. Initial conditions. The geometry of the striped magnetic field is shown in Fig. a). The longi- 
tudinal Px and transversal py component of the momentum of one species are seen in Fig. b) and c) 
respectively. The temperature and the drift speed for one species (the speed of the other species is 
equal and opposite) are shown respectively in Fig. d) and e). 



For the results presented in this section, the simulation box is divided into 2x10^ cells and we 
used 8 X 10^ particles for each species, i.e. 8 x 10^ electrons and 8 x 10^ positrons. Therefore, on 
average, there are only 4+4 particles per cell. Although this number seems rather small, we checked 
that increasing the number of particles per cell do not improve the accuracy of our results. Thus, 
we performed the whole set of simulations with this resolution. To support our statement, we will 
show examples with 40+40 particles per cell. The particle number density in the hot unmagnetised 
part, i.e. in the sheets, is five times higher than in the cold magnetised part, Whi = 5 nd. The time 
step is chosen such that Af - 0.5. The resulting upstream Larmor radius, deduced from this 
value and from Fi - 20 according to Eq. ( l28T l, is equal to ro = 40 cells. Half a wavelength of 
the striped wind is /i = 2500 cells and the relative current sheet thickness is 6i =0.1. At the left 
boundary of the simulation box, no incoming plasma is injected. The upstream plasma is simply 
flowing to the right and leaves an empty space behind it with a constant magnetic field advected 
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at the flow bulk speed Fi . In order to keep meaningful results, the simulation has to be stopped at 
the time ff - L/2c before the electromagnetic wave starting at f = at the right wall reaches this 
vacuum region, propagating from left to right starting from x = at f = 0. 

In all runs, the simulation box starts with a inhomogeneous plasma made of electrons and 
positrons following a 2D relativistic Maxwellian distribution function in their proper frame given 
by 



with relativistic temperature ©j = kg TJiris and Lorentz factor 7=^1+ {pi + Py)/ml c^. The 
temperature 0^ is obtained from the pressure equilibrium condition Eq. ^ and shown in Fig.[3ji)- 
A^j is the particle number density in the proper frame, m, the mass of a particle of each species. 
Because they are equal for electrons and positrons, we note = m. We have to distinguish between 
3 reference frames 

- R^: the species frame or proper frame of reference for each species, electrons and positrons, in 
which the distribution function is given by the 2D relativistic Maxwellian Eq. i5% ; 

- Rv,: the wind frame of reference in which the leptons are counter-streaming with velocity + Us ; 

- ^0- the observer frame or simulation box frame or lab frame of reference in which the wind 
(the current sheets) is propagating in the positive x-direction with a Lorentz factor Fi . 

To upload the initial distribution function of each species in the lab frame, we have to make two 
successive Lorentz transformations, the first one leading from R^ to /?« and the second one leading 
from /?„ to Rq. 

We performed many runs with different initial Lorentz factor of the wind and different mag- 
netisation. We give a typical sample of runs demonstrating fuU dissipation of the magnetic field or 
no dissipation at all. 

Because the velocity space is only two-dimensional in our simulations, see the distribution 
function Eq. ( |52] l, the adiabatic index for the relativistic plasma is y = 3/2 instead of the usual 
4/3. This different index will affect the parameters downstream the shock like the temperature, the 
Lorentz factor and the magnetisation. The general jump conditions for arbitrary adiabatic index are 
derived in appendixlAl These formulae will help us to make quantitative comparisons between the 
analytical results presented before (and made for y = 4/3) and the PIC simulations. 

4.1. Negligible dissipation 

Here we show details of a typical example where magnetic dissipation is negligibly small. The 
magnetisation between the sheets is set to cr = 3. It corresponds to the definition of crj in Eq. ([T]). 
The true magnetisation Eq. (|53] |. equal to cr 2.7, differs from the one in the cold part cr because 
in our PIC simulations, we need to resolve the structure of the sheet (sheet thickness larger than the 
Larmor radius) therefore the contribution from the cold unmagnetised part is not negligible. The 
true average magnetisation is actually roughly 10% lower than the one in the cold part due to the 
relative current sheets thickness 6i =0.1. The true magnetisation, averaged over one wavelength 
of the striped wind , is given in the rest frame of the downstream plasma by 



e 



iy-i)/@. 



MX,p^,Py) = — 



2n m?c2 0,(l H-0,) 



(52) 
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For a fluid with adiabatic index y = 3/2, the average enthalpy density (in the rest frame) is given 
by 

(enthalphy density) = 3/2 (internal energy density) - 1 /2 (particle number density) (54) 

The internal energy includes the rest mass energy as well as the kinetic energy of the particles. 
This formula is only valid in the rest frame of the fluid. In our computations, the frame of the 
simulation box corresponds to the rest frame of the downstream plasma, therefore Eq. ( l54l l applies 
to this plasma but not to the upstream one which moves with a Lorentz factor Ti with respect 
to the simulation box. Thus, we have first to transform back to the rest frame of the upstream 
plasma and then apply Eq. (|53l l and (|54] i. The upstream rest frame is easily found because we 
impose the Lorentz factor of the upstream plasma with respect to the simulation box which is also 
the rest frame of the downstream plasma. Therefore, the Lorentz transform in the upstream frame 
involves simply Fi . Note however that this is correct only well upstream and not in the shock 
itself. The bulk Lorentz factor diff'ers from Fi whenever a precursor arrives. Nevertheless, we are 
only interested in the average quantities when the front shock has passed through the plasma. The 
precise values within the shock front are not significant for our study. Note that the true local 
wavelength upstream as well as downstream the shock, necessary for the averaging <> in Eq. (l53T l 
and (|54] |. is determined by looking for the locations where the magnetic field vanishes. Indeed, 
due to compression of the plasma after shock crossing, the length of the stripes downstream are 
reduced by a factor y6i/yS2(ji:) > 1 compared to their value upstream, according to Eq. (flST l. Note 
also that the downstream wavelength can vary slightly (a few percent) from one period to the next 
one, due to some perturbations. These perturbations propagate also in the upstream flow, causing 
some small changes in the upstream wavelength too. All these variations are taken into account 
when computing the averaged quantities presented in the figures. 

The results of the run are summarised in Fig.|4]for a final time of simulation ff - L/2 c - 5x10^. 
This corresponds to the time needed by an electromagnetic wave to propagate from the wall at the 
right boundary, starting at f = and arriving at the middle of the simulation box L/2 at the time ff . 
The shock front is very sharp with a very small thickness and located at approximately x/L ^ 0.62. 
It propagates in the negative x direction, starting from the wall located at the right boundary as 
before. Looking at the structure of the magnetic field. Fig. |4^), the stripes downstream are per- 
fectly preserved. They are still clearly recognisable in Fig. ^) and c). They are only subject to a 
compression of a factor roughly equal to two. Thus the wavelength of the wind in the downstream 
frame has been divided by two whereas the magnetic field strength, due to magnetic flux conserva- 
tion, increases by a factor two. The mean particle density in the shocked plasma. Fig. 15), reaches 
a value close to two, in accordance with the magnetic field compression ratio of two. Because we 
average over one wavelength of the striped wind, the diff'erence in hot and cold density is smoothed 
out. 

The factor two can be explained as follows. The shock propagates with velocity close to c. 
In the shock frame, the magnetic field strength and the particle number density remain nearly the 
same, B2 = Bi and n2 = "i (in the high cr Umit the shock is weak!). However, in the simulation 
frame, which is the downstream frame, they vary because of Lorentz transformations. Indeed, let's 
note n 1 /2 the density of the upstream plasma, proper density n j , as measured in the rest frame of the 



18 Petri and Lyubarsky: Magnetic reconnection in a striped pulsar wind 

downstream plasma. Let also ri/2 be the relative Lorentz factor of the incoming flow as measured 
in the shocked plasma. It is easy to show that a Lorentz transformation gives 

ri/2 = rir2(i-/3iy82) (55) 

Thus the upstream density as measured in the downstream plasma is 

ni/2 = ri/2n'i (56) 
With help of the jump condition for particle number density conservation, we found 

^ - (57) 



ni/2 ri/2r2^2 

In the ultra-relativistic limit, ySi — > 1, we get 

«2 i+yS2 



(58) 



«l/2 Pi 

Therefore the factor 2 when y62 ~ 1. Actually, because the magnetisation is not very high, cr\ - 
3, in the shock frame, the downstream plasma propagates with a Lorentz factor r2 ~ V5 o"i/4 
(see appendix IaIi. corresponding to a speed y62 = 0.86. Inserting this value in Eq. (ISSl l. we found 
«2/"i/2 ~ 2.15, very close to the average density in the shocked plasma, see Fig.|4^). Therefore, the 
shock front propagates in the negative jc-direction with the speed ^62. Starting from f = at x = L, 
at the end of our simulation, it traveled a distance Ax = j02 c/2 * ff 0.43 L. Therefore the shock 
front should be located at xjL - 0.57. However, the magnetisation at the initial time close to the 
wall is very weak, close to unity. Therefore, the shock speed is smaller at earlier times. That is why 
the shock front is located at only x/L « 0.62 as can be checked by inspecting Fig.|4^) and b) and 
c) where a sharp transition is observed between the upstream and downstream flow. 

Another way to estimate the speed of the shock is to compare the position of the shock front at 
the final snapshot, (x/L a; 0.62) and at the penultimate snapshot (x/L a; 0.664). Moreover, the time 
between these two snapshots is Af - LjlQc. Therefore, the speed of the shock front is roughly 
j82 ~ 0.044 L/(L/20) ~ 0.88, close to flie estimate made above. 

The particles downstream are not thermalised. Indeed, the particle momentum which is mainly 
directed along the propagation axis of the wind in the upstream flow, py <K px, is only weakly 
disturbed by the shock. The magnitude of is not altered by the shock. Nevertheless, after crossing 
the shock, both longitudinal and transversal momenta are of the same order of magnitude, Py ~ Px- 
The mean particle Lorentz factor, y, is roughly increased by a factor of two. Fig. |5] Although 
the transition for the transverse momentum component is very sharp, the stripes are still clearly 
identifiable, Fig.l^^) and c). The magnetisation downstream is increased by a factor roughly equal 
to two, as expected from the analysis of Sect. [3] Fig. H)!), demonstrating that the magnetic field 
does not dissipate. Contrary to the run presented in the next section, there is no electromagnetic 
precursor, the flow is not perturbed by the shock front. The situation is very similar to the ideal 
relativistic MHD shock. The magnetic field lines, frozen into the plasma, have to follow the motion 
imposed by the matter. The shocked plasma is compressed, and due to the frozen magnetic flux, 
the stripes are also compressed in the same ratio. 

Increasing the number of particles per cell will not improve the accuracy. Indeed, we performed 
another run with ten times more particles per cell (40 + 40), see Fig. |6] and compare with Fig. ID 
Therefore, using only 4 + 4 particles per cell is justified. 

Finally, we show an example with zero average magnetic field a - Oin Fig. [T] to prove that this 
parameter a does not affect much the flow downstream provided a 1. 
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Fig. 4. Example of run showing no magnetic reconnection at the termination shock. The Lorentz 
factor of the wind is T = 20 and the magnetisation cr = 3. The magnetic field is shown in Fig. a), 
the particle longitudinal and transverse momentum component in Fig. b) and c) respectively, the 
intermediate, t - 2.5 xlO^ (dotted line) and final, t - 5x10^ (solid line) downstream magnetisation 
in Fig. d) and the intermediate (dotted line) and final (solid line) average particle density number 
downstream in Fig. e). All quantities are plotted against the normalised abscissa x/L. 




Fig. 5. Mean particle Lorentz factor averaged over one period of the striped wind at an intermediate, 
J = 2.5 X 10"* (dotted line) and final, / = 5 x 10"* (solid line) time of simulation. 

4.2. Full dissipation 

In this paragraph, we discuss in detail a typical example of full magnetic reconnection in the striped 
wind. The magnetisation in the cold magnetised part, i.e. between the sheets, is set to cr = 45. The 
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Fig. 6. Same as Fig.|4]but with ten times more particles per cell, i.e. 40 positrons and 40 electrons 
per cell. The results are nearly the same as in Fig.|4l However, we stopped the run when made sure 
that everything proceeds like the run with 4 + 4 pairs so that the position of the shock is different. 



true average magnetisation, ctj ^ 40.0, is here again roughly 10% lower than the one in the cold 
part due to the relative current sheets thickness 6i =0.1. 

The only change we made compared to the previous case is to take a higher magnetisation. This 
has dramatic consequences on the flow downstream as we will see now. 

The simulation results are shown in Fig. [8] for the final time tf = L/2c = 5 x 10"*. At this 
snapshot, the shock front is located at approximately x/L « 0.75 (see below how this location 
was determined). It propagates in the negative x direction, starting from the wall located at the 
right boundary x - L. Looking at the structure of the magnetic field. Fig. [8^), it is clearly seen 
that the stripes, after crossing the shock are destroyed. Nevertheless, because the average magnetic 
field does not vanish, a small DC component remains, of the order of a = 0.12. The particles down- 
stream are fully thermalised. Indeed, in the upstream flow corresponding to x/L < 0.78, the particle 
momentum is mainly directed along the propagation axis of the wind, longitudinal momentum 
depicted in Fig. |8j)), with non negligible transverse momentum py only within the current sheets. 
Fig- lit)- The upstream longitudinal component p^ is by orders of magnitude larger than the trans- 
verse component py, py <k p^. The particle momentum distribution function is randomised when 
crossing the shock discontinuity, i.e. in the region defined hy x/L > 0.75, and therefore downstream 
the two component of the momentum are of the same order of magnitude, py x p^ and by orders 
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Fig. 7. Same as Fig.|4]but with zero average magnetic field, a - Q. The results are nearly the same 
as in Fig.m 



of magnitude larger than before the shock. The enhancement in the momentum of the particle in 
both directions, /?x and py, is explained by the Poynting flux dissipation. Magnetic energy is con- 
verted into particle kinetic energy. There are no more spikes in the px graph proving that the stripes 
disappeared. The magnetisation downstream is almost zero, solid curve in Fig.|8}l), indicating that 
the alternating magnetic field has completely dissipated into particle heating as expected. In that 
case, the flow downstream is hydrodynamical and the compression ratio close to three. The mean 
particle density in the shocked plasma, solid curve in Fig.|8^), indeed reaches a value close to three. 

The factor three can be explained in the following way. Assuming full dissipation of the mag- 
netic field, the upstream flow is strongly magnetised whereas the downstream plasma is purely 
hydrodynamical. For an ultra-relativistic gas, the adiabatic index is y = 4/3. This is true for parti- 
cles evolving in a three dimensional velocity space in which they possess three degree of freedom 
in translation motion. However, in our simulations, the distribution function, Eq. (|52T i. is only two- 
dimensional in velocity space and the corresponding adiabatic index is thus y = 3/2. Solving the 
MHD jump conditions for this 2D plasma, it is easily found that the downstream plasma velocity 
in the rest frame of the shock is /?2 = 1 /2 (2D velocity space) instead of the traditional /32 - 1/3 
(3D velocity space). Using Eq. (|58T l, the ratio of the particle number density expressed in the down- 
stream frame (simulation box frame) is /n 1/2 - 3 as expected. Again, the front shock propagates 
in the negative x-direction with the speed ySa = 1/2. Starting from f = at x = L, at the end of our 
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simulation, it traveled a distance Ax - c/2 tf - L/4. Therefore the shock front should be located at 
x/L - 0.75. However, because of some small retardation effect to initiate the shock front, it is seen 
to lie at x/L ^ 0.78 by inspecting Fig. [8^) and b) and c). 

It is also interesting to note that the incoming flow is already perturbed before entering the 
shock. The magnetisation gradually decreases and becomes already very small in the region just 
before the shock discontinuity. The stripes start to be significantly disturbed akeady at x/L x 0.6. 
This phenomenon is probably due to an electromagnetic precursor, i.e. and electro magnetic wave 



generated at the shock fron t and propagating downwards into the upstream plasma dHoshino et al 



1992 : 



Gallant et al 



19921). Moreover, from the momentum component p^, we conclude that the 



flow is heated already well upstream of the shock. If the strong electromagnetic precursor would 
be generated immediately at the initial time f = 0, we would expected it to influence half of the 
simulation box from x/L = 0.5 Xo x/L - 1 at the final time ff. However, the amplitude of the 
precursor is initially small because due to the chosen initial distribution of the flow parameters 
(FigO, the energy of the plasma entering the shock is initially small. The strong shock is formed 
after some time (about 0.05L/c) when the highly relativistic plasma enters the shock. Thus a strong 
enough precursor is generated not from the beginning of the simulations and therefore the flow is 
significantly perturbed at x/L > 0.55. The effect of a precursor was not included in our analytic 




Fig. 8. Example of run showing the magnetic reconnection at the termination shock. The Lorentz 
factor of the wind is F = 20 and the magnetisation cr - 45. The same quantities as in Fig. |4] are 
plotted. 
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analysis. It is performed with the assumption that the dissipation is weak; in this case, there is no 
significant precursor In any case, our analysis is valid provided the downstream flow is settled into 
a striped wind. Then we apply the conservation laws to flows well upstream and well downstream 
of the shock and processes within the "black box" containing the shock and possible precursors do 
not affect the results. 

In case of high upstream magnetisation, the flow in the wind is dominated by the dynamics 
of the magnetic field, i.e. gaseous pressure negligible compared to magnetic pressure. Therefore 
the shocked plasma has to reorganise in such a way to keep the striped structure. If this plasma is 
unable to maintain the current required by the magnetic field, the latter will dissipate. 

Finally, here again, we show an example with zero average magnetic field ff = in Fig.|9] One 
sees that the Poynting flux dissipates completely in this case also. However there is no sharp shock 
transition in this case. When the average magnetic field is nonzero, the shock is mediated by the 
Larmor rotation of particles and the width of the shock transition is roughly equal to the Larmor 
radius in the average field, ro/a. In our simulations this is 4Q/a cells so that even at a - 0.1 the 
shock transition is narrow. When a = 0, the particle free path is determined only by scattering on 
magnetic fluctuations therefore the shock transition could be rather wide. 




Wife 



Fig. 9. Same as Fig.|8]but with zero average magnetic field, a = 0. 
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4.3. Empirical law 

In the last two paragraphs, we gave two examples of runs, the former demonstrating that the striped 
wind can be preserved when some conditions are fulfilled and the latter showing a case of full 
magnetic reconnection. We performed many simulations with different parameters, changing the 
magnetisation, the Lorentz factor, the Larmor radius in order to segregate between the two regimes. 

The results of the full set of simulations can be summarised in plots showing the ratio of the true 
downstream to the true upstream magnetisation 0-2! cr\ as well as 0-2, Fig. [10] To sum up, we obtain 
the curves represented in Fig. [10] where these quantities are plotted against the parameter 77 - 
II ro cr'l^ and 77 ^[a■ for fig. a) and b) respectively and for different initial magnetisations, as was 
done in the semi-analytical results presented in Fig.|2]e) and f). 
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Fig. 10. Empirical law for magnetic reconnection in the striped wind obtained by ID PIC sim- 
ulations. The ratio of the downstream to the upstream true magnetisation crtja-y is plotted 
against Z/ro cr3/2 in fig a) and the downstream magnetisation against Ijrocr in fig b) for differ- 
ent values of the theoretical upstream magnetisation cr from mildly cr = 1 to highly magnetised 
flow cr = 75. 
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Fig. 11. Same figures as Fig[2t) and f) respectively fig. a) and b) here, but for an adiabatic index 
y - 3/2 instead of 7 = 4/3. 

For the largest values of 77, no reconnection is observed as we would expect from the analytical 
study discussed in Sect. [3] In the no dissipation limit, 77 » 1, the ratio crij(T\ reaches values close 
to 1.7, in accordance with the expected 5/3-ratio for the ideal MHD shock with 7 = 3/2 (see 
appendix lAb. As this limit is achieved at the same 77 for any cr, we retrieve the condition Eq. ( l42b . 
According to Fig[TO]b) all the presented curves, with the exception of the one corresponding to 
(J - 1, go to zero at the same condition 

— < 3. (59) 
roCT 

This confirms the analytical criterion|45]found at the condition cr » 1. 

To allow detailed comparisons of the PIC simulations with the theoretical model presented in 
sections|2]and[3] we plotted in Fig[TT]the quantities 0-2/0-1 and 0-2 calculated from the theoretical 
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model with the adiabatic index 7 = 3/2 (the same quantities were plotted in Figs|2^) and f) 
for J - 4/3). Theoretical jump conditions for arbitrary y are presented in Appendix B. Comparing 
Fig.[TT]a) and Fig. [Tola) on one hand, andFig.fTTIb) andFig.fTOlb) on the other hand, the parameter ^ 
(the ratio of the sheet width and the particle Larmor radius), introduced at the end of Sect. |2l 
can be easily estimated. Indeed, looking at figfTOb). we find that no dissipation occurs whenever 
//n, 0-3^2 ^ whereas figEJi) gives llrQ^fr"'^ * 10 so that ^ * 10. On the other hand, figHOb) 
shows full dissipation at l/rocr ^ 3. Comparing with fig. [TTb) which gives l/ro^cr ^ 0.5, we get 
^ ss 6. This means that ^ a: 6 - 10. 

Finally substituting Eqs.([T]i and ( |27T i into Eq.(|59]l, the parameter region for full dissipation is 
found as 

y"oc|?|^<3 (60) 

Note also that the condition for full dissipation Eq. ( |60] | is independent of the upstream Lorentz 
factor F]. 



5. Application to pulsar winds 

We apply Eq. ( |60l l to pulsar winds. The idea is to find the limiting radius beyond which the alter- 
nating field dissipates at the termination shock. Because all quantities refer to the upstream plasma, 
we drop the subscript 1 . 

Beyond the light cylinder, we assume that the outflow is radial and propagating in the radial 
direction. The magnetic field is predominantly toroidal and decreases with radius as 
Rl 

where R]^ - cjQ. is the radius of the light cylinder, Q. the angular velocity of the neutron star and 
Bl the magnetic field strength at the light cylinder. The particle number density falls off as 

« = "L ( (62) 

where the density at the light cylinder hl is conventionally expressed via the so called multiplicity 
factor K such that 

hl = 2 £() K — — (63) 

\q\ 

Half of a wavelength in the striped wind is given by I - nR]^. Substituting these relations into 
Eq.(l60ll. one finds that full magnetic dissipation occurs at the termination shock whenever the 
shock arises at the distance larger than 

We suppose that the termination shock is stationary in the observer (pulsar) frame. 

The value of the multi plicity, k, is rather uncertain. The available theoretical models 

ousands. The observed synchrotron emis- 



jHibschman & Arons , 



sion from the Crab nebula places wide limits on k, from k ~ 10^, at the assumption that only high 
energy electrons (those emitting the optical and harder radiation) are injected now into the nebula, 
io K ~ 10^ if the radio emitting electrons are also injected now but not only in the early stage of 
the pulsar history. In any case the condition Eq. ( |64] l is easily satisfied for plerions because the 
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termination shock is located at radii 10'"^ - 10'* m ~ 10** - 10"^7?l from the pulsar For instance, 
in the Crab nebula, the termination shock is located at /?shock/^L = 3 x 10^ » k. Thus for these 
pulsar wind nebulae, the magnetic energy is easily dissipated at the termination shock so that the 
pulsar wind could remain highly magnetised even in the upstream region close to the shock front. 
All the available observation limits on cr are obtained from the analysis of the plasma flow and 
radiation beyond the termination shock, the necessary upstrea m cr being calculat e d from the ideal 



MHD j ump condit i ons, s ee for instance the recent papers by 
feooi; 



Torii et al 



Petreetal 



mm -. 



Safi-Harb et al 



(120021) ■ Our results show that even if the pulsar wind remains Poynting domi- 
nated until the termination shock, magnetic dissipation at the shock front would suffice to efficiently 
convert the electromagnetic energy to the plasma energy. This point of view is an alternative to the 
cr-problem. 

In binary pulsars, interaction of the pulsar wind with the companion star results in forma- 
tion of an intra-binary s hock very close t o the neutron star, only tho usands of light cylinder radii. 



Arons & Tavani 



( ll993l) : lTavani & Aronsl(ll997h : ICheng et alj ( 120061) claim fliat flie X-ray emission 
observed from the binary pulsars PSR 1957+20 and PSR 1259-63 is explained if the magnetisation 
downstream of the shock is small so that the Poynting flux is already converted into the plasma 
energy. Then we could estimate from Eq.(l64li the upper limit of the pair production factor k for 
these pulsars. 

Indeed, first consider the pulsar PSR 1957+20. It has a rotation period of P = 1.607 ms, corre- 
sponding to a light cylinder radius of = 76.7 km. The intra-binary shock arises at a distance less 
than the orbital separation to the companion star, so that we get i^shock S 1.7 x 10* km. According 
to Eq. ( l64b . in order to expect a significant dissipation of the magnetic field at the shock, the mul- 
tiplicity factor should then satisfy k < {I In) /?shock/^L ~ xlO**. 

Next consider the pulsar PSR 1259-63. It possesses a period of P = 47.7 ms, corresponding 
to a light cylinder radius Pl - 2279 km. Because here again the shock should arise at a distance 
smaller than the orbital separation, we found Pshock ^ 3.9 x 10** km. Magnetic reconnection implies 

K<%X 

Another inter e sting case is the famous double pulsar PSR J0737-3039. According to 



McLaughUn et al.l (120041) . the radio emission from the pulsar B is modulated with the period of 



pulsar A. This strongly supports the idea that the magnetosphere of pulsar B is disturbed by the 
striped wind emanating from pulsar A. The condition that the striped structure has not been erased 
at the distance separating both pulsars, imposes a lower limit on the multiplicity factor k in this 
case. Let us estimate this limit. The period of pulsar A is - 22.699 ms corresponding to a 
light cylinde r radius of Pl - 1083 k m. The separation between the two pulsars is approximately 



700, 000 km ( Lvne & Kramei 



2OO5I) . Then die lower Hmit for the ratio in Eq. is P/Pl ~ 646. 



We conclude that the lower Hmit for the pair production multiplicity in pulsar A is a: > 310. 



6. CONCLUSION 

In this paper, we studied the magnetic reconnection process at the termination shock in the pul- 
sar striped wind. Using the jump conditions for the averaged quantities (over a wavelength) in 
the striped wind, we derived a simple analytical criterion for magnetic dissipation, see Table [1] 
However, we had to introduce a free parameter ^, Eq. ( [38l ), entering the dissipation condition. It 
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relates the downstream thickness A2 of a current sheet to the downstream Larmor radius ro by the 
following expression A2 = ^ tq. ^ is the only free parameter in our model. In order to estimate this 
parameter, we performed ID relativistic PIC simulations and found that ^ ^ 10. 

Knowing the important flow parameters in the incoming plasma which are half the wavelength I, 
the magnetisation cr and the Larmor radius ro, we are able to predict the percentage of magnetic 
reconnection in the striped wind when crossing the termination shock, according to the following 
law 

- for l/rocr < 3, full dissipation occurs, the flow downstream is purely hydrodynamical with 
Lorentz factor close to unity Fa ~ 1, particle are thermaUsed and heated to relativistic tempera- 
ture. The magnetic field has completely disappeared, except for the DC component (i.e. average 
magnetic field) ; 

- for cr < (1/12 ro)^^^, no magnetic reconnection exists. The striped wind structure is preserved 
downstream, it is just compressed. The downstream Lorentz factor is the same as for ideal 
MHD, Tz = Vo" ; 

- 3 < l/rocr < 12 Vo", the magnetic field is partially dissipated. The stripes are weakened and 
the flow is decelerated to a downstream Lorentz factor Fz = 1/(12 ro cr). 

We applied the condition for full magnetic dissipation, Eq. ( |64] |. to pulsar wind nebula and 
binary pulsars. Because in plerions the termination shock is located at radii lO''* - lO'^m ~ 
10** - 10'" /?L from the pulsar, magnetic reconnection is easily achieved. This conclusion could 
resolve a long-standing difficulty with transformation of the electro-magnetic energy of the pulsar 
wind into the plasma energy in the nebula (th e so called cr- problem). First of aU, both observa- 



tions of the X-ray tori jWeisskopf et all 120001) and theoretical models dBogovalovi 119991) suggest 
that most of the energy in the pulsar wind is transferred in the equatorial belt where the mag- 
netic field is predominantly alternating. We see now that even though dissipation of the alternat- 



ing fields in the wind is h ampered by relativistic slowing-down of time dLvubarskv & Kirk , 



2001 



Kirk & Skiaeraasenll2003h . the electro-magnetic energy is readily released at the terminating shock. 

Our model could be applied also to pulsars in binary systems, where interaction of the pulsar 
wind with the companion star results in a formation of a shock relatively close to the pulsar. In this 
case, the presence or absence of magnetic dissipation imposes a strong constraint on the pair mul- 
tiplicity factor K. For pulsars PSR 1259-63 and PSR 1957H-20, where full dissipation is expected, 
the upper limit is roughly k < few x 10"*. For the binary pulsar PSR 0737-3039, we do not expect 
full reconnection and therefore the lower limit is found to be a- > 310. 

The generation of electromagnetic waves at the shock front and their propagation in the up- 
stream plasma will affect the flow already before entering the discontinuity. This aspect of the 
reconnection in the termination shock is left for future work. 
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Appendix A: Jump conditions for arbitrary adiabatic index y 

A.1. Ideal MHD 

In this appendix, w e give the general expressions for the ideal MHD jump conditions for arbitrary 



adiabatic index y, dGallant et al 



1992; 



Hoshino et al 



19921) . Let's assume that the equation of 



state for the ultrarelativistic plasma is given by 
P 



(A.1) 
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for an arbitrary index j. The enthalpy is therefore given by 
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Kennel & Coroniti 



IP 

w = 

7-1 

Following the procedure of 
for the 4-velocity gives 



2 _ (4-y)ycr]+4(l+(j-lf)cri+4(y-lf 



(A.2) 



(Il984al) . the solutions to the MHD jump conditions 



(A.3) 



^16 r(r - 2)3 cri (cTi + 1) + [(r - 2)2 crl-4(cri + l) (en + (y - l)^)]^ 
- 8y(2-y)(cri + l) 

Starting from the 4-velocity, we can compute the Lorentz factor, the downstream temperature and 
the magnetisation by 



kBT2 ^ y-l 1 

r ^2 



y-l n2 Ti m c 

(T2 = O"! 

y ni 12 kBT2 
For 7 = 4/3 we retrieve the usual result 



8 cr2 + 10 0-1 + 1 ± ^64cr| + 128cr3 + 84cr2 + 20cri + 1 
2 16(cri + l) 

For arbitrary y, the expansion in l/cri to the second order leads to 



y-4 



"2 4(7-2)^^' 4(72-67 + 8) 
The corresponding Lorentz factor is to the same order 

J, ^1 l (cri+3)72-8((ri+2)7 + 8(2cri+3) 
^ 2 ^| 72-67 + 8 

This leads to a downstream temperature 



7" - 8 7 ■ 



fcfiTa _ 7 - 1 y- 2 1 
T\mc^ 4-7^7-4 ^J(^l 

a density jump 

^ = 1+22:^^ 

«! 7 - 4 CTi 

and a magnetisation 



/4-7 ^ (y-2)2 . 

0-2 = 2 — — 0-1 



2- 



(3r'-87 + 8) 1 \ 



(7 - 4)2 7^ CT] 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



(A.8) 



(A.9) 



(A. 10) 



(A. 11) 



(A. 12) 



(A.13) 



\ r r(r-4)cri, 

InTable pV.ll we give the numerical values for the different coefficients for 7 = 4/3 and 7 = 3/2. 
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Quantity 7 = 4/3 7 = 3/2 

n2/ni 1 + 1 + 



2^1 

0-2/0-1 2 + f + 



Table A.l. Downstream parameters for the two adiabatic indexes y - A 13, 3 12, expanded to first 
order in the upstream magnetisation 1 /tri . 



A.2. Striped wind 

We now give the general formulae for the jump conditions in the striped wind for arbitrary adiabatic 
index y. Following the same procedure as in Sect |2] using the general expression described in the 
previous section, the average conservation of particle and energy-momentum read. 



j6i [(I -61) Hci + ^1 "hi] = jSi [(1 - 62) nc2 + S2 rihi] 



P\ tic 



1-61 + 



2-r 



2(r-i) 



6i + l\cr 



■y 



2(r-i) 



t 2 



r 



+ 11 + J—^ 62 



2(y- 1) / Hci 



nc2 



■CTl 



r 



r- 1 1 _^ 



4 - y 

nc2 0-i 



+ 1 + 



■7 



2(7-1) 



(A. 14) 
(A. 15) 



4 - 7 2 «ci 
The prescription for the current sheet thickness downstream is 

ySl Hcl '"0 



^2 I — 0-1 + 

"cI / 



(A.16) 



52 =^ 



7-1 _^ nc2 



4-7 2«ci 



r2y82nh2 1 1 



(A.17) 



